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Spatial entanglement using a quantum walk on a many-body system 

Sandccp K. GoyaQ 
The Institute of Mathematical Sciences, CIT campus, Chennai 600 113, India 

C. M. ChandrashekaiQ 

Institute for Quantum Computing, University of Waterloo, Ontario N2L 3G1, Canada and 
Perimeter Institute for Theoretical Physics, Waterloo, ON, N2L 2Y5, Canada 

The evolution of a many-particle system on a one-dimensional lattice, subjected to a quantum 
walk can cause spatial entanglement in the lattice position, which can be exploited for quantum 
information/communication purposes. We demonstrate the evolution of spatial entanglement and 
its dependence on the quantum coin operation parameters, the number of particles present in the 
lattice and the number of steps of the quantum walk on the system. Thus, spatial entanglement 
can be controlled and optimized using a many-particle discrete-time quantum walk. 



I. INTRODUCTION 



Entanglement in a quantum state has been the fundamental resource in many quantum information and com- 
putation protocols, such as cryptography, communication, teleportation and algorithms [H, 0]. To implement these 
■ protocols, generating an entangled state is very important. Similarly, studies on the interface between condensed 
matter systems and quantum information have shown entanglement as a signature of quantum phase transition [H-Q ■ 
4—> ■ To understand the phases and dynamics in many-body systems an analysis of entanglement in many-body systems is 
very important. Hence, various schemes have been proposed for entanglement generation in quantum systems 
and for understanding entanglement in many-body systems (Toj . Quantum walk (QW) is one such process in which 
Q-i an uncorrelated state can evolve to an entangled state and be used to analyze the evolution of entanglement [ill, EH ■ 
i— —i \ The QW, which was developed as a quantum analog of the classical random walk (CRW), evolves a particle into an 
entanglement between its internal and position degrees of freedom. It has played a significant role in the development of 
quantum algorithms (l3j . Furthermore, the QW has been used to demonstrate coherent quantum control over atoms, 
^ j i quantum phase transition 14 1, to explain the phenomena such as breakdown of an electric field-driven system (l5j 
' and direct experimental evidence for wavelike energy transfer within photosynthetic systems [Ttl Il7j . Experimental 
\Q , implementation of the QW has also been reported fl8l - [2l| . and various other schemes have been proposed for its 
physical realization [22l426j . Therefore, studying entanglement during the QW process will be useful from a quantum 
* , information theory perspective and also contribute to further investigation of the practical applications of the QW. 
In this direction, evolution of entanglement between single particle and position with time (number of steps of the 
0"\ \ discrete-time QW) has been reported [l2j . 

In this paper, we consider a multipartite quantum walk on a onc-dimcnsional lattice and study the evolution of 
spatial entanglement, entanglement between different lattice points. All the particles considered in the system are 
identical and indistinguishable with two internal states (sides of the quantum coin). Spatial entanglement generated 
using a QW can be controlled by tuning different parameters, such as parameters in the quantum coin operation, 
number of particles in the system and evolution time (number of steps). To quantify entanglement in the system we 
are using Meyer- Wallach multipartite entanglement measure. 

In Sec. UH we describe single-particle and many-particle discrete-time QWs. In Sec. IIII1 entanglement between a 
particle and position space and spatial entanglement using single- and many-particle QWs arc discussed. In Sec. IIV1 
we present the measure for spatial entanglement of the system using the Meyer- Wallach global entanglement measure 
scheme for particles in a one-dimensional lattice and in a closed chain (n— cycle). We also demonstrate control over 
spatial entanglement by exploiting the dynamical properties of the QW. We conclude with the summary in Sec. fVl 
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II. QUANTUM WALK 

Classical random walk (CRW) describes the dynamics of a particle in position space with a certain probability. The 
QW is the quantum analog of CRW-developed exploiting features of quantum mechanics such as superposition and 
interference of quantum amplitudes [27l - l29j . The QW, which involves superposition of states, moves simultaneously 
exploring multiple possible paths with the amplitudes corresponding to the different paths interfering. This makes 
the variance of the QW on a line to grow quadratically with the number of steps which is in sharp contrast to the 
linear growth for the CRW. 

The study of QWs has been largely divided into two standard variants: discrete-time QW (DTQW) (29l - [3l| and 
a continuous-time QW (CTQW) [32|. In the CTQW, the walk is defined directly on the position Hilbert space H p , 
whereas for the DTQW it is necessary to introduce an additional coin Hilbert space H c , a quantum coin operation 
to define the direction in which the particle amplitude has to evolve. The connection between these two variants and 
the generic version of the QW has been studied [HI, Hi[. However, the coin degree of freedom in the DTQW is an 
advantage over the CTQW as it allows control of dynamics of the QW [H, [36|. Therefore, we take full advantage of 
the coin degree of freedom in this work and study the DTQW on a many-particle system. 



A. Single-particle quantum walk 

The DTQW is defined on the Hilbert space H. = H c ® H p . In one dimension, the coin Hilbert space H c , spanned 
by the basis state |0) and |1), represents two sides of the quantum coin, and the position Hilbert space Ji p , spanned 
by the basis states \i/jj), j € Z, represent the positions in the lattice. To implement the DTQW, we will consider a 
three-parameter U(2) operator C^ } g^ of the form 

= / e^cos(e) e z <sm(6) \ , , 

°f AC ~ I S in(6») - e -« cos(0) ' 



e 



as the quantum coin operation [36| . The quantum coin operation is applied on the particle state (Cf e,£ ® 1) when 
the initial state of the complete system is 

|*in) = [cos(*)|0) + e ir > Bin(<5)|l>] ® |Vo>- (2) 

The state cos(<$)|0) +e 171 sin(5) 1 1) is the state of the particle and |^o) is the state of the position at the lattice position 

j = o. 

The quantum coin operation on the particle is followed by the conditional unitary shift operation S which acts on 
the complete Hilbert space of the system: 

S = exp(-io- z <g> PI), (3) 

where P is the momentum operator, o~ z is the Pauli spin operator in the z direction and / is the length of each step. 
The eigenstatcs of a z are denoted by |0) and |1). Therefore, S, which delocalizes the wave packet over the positions 
(J — 1) and (j + 1), can also be written as 

s = |o)(o| ® J2 \^-i)(i>j\ + |i>(i| ® E (4) 

The process of 

W L 8£ = S(C 6A( (g) 1) (5) 

is iterated without resorting to intermediate measurement to help realize a large number of steps of the QW. The 
parameters S and 77 in Eq. ([2]) can be varied to obtain different initial states of the particle. The three parameters 
£, 8 and £ of C^^x can be varied to choose the quantum coin operation. By varying parameter 9 the variance can 
be increased or decreased according to the functional form, a 2 w (1 — sin(6>))£ 2 , where t is the number of steps of the 
QW, as shown in Fig. [TJ 

Biased coin operation and biased QW: The most widely studied form of the DTQW is the walk using the Hadamard 
operation 
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FIG. 1: (color online) Spread of the probability distribution for different values of 8 using the quantum coin operator Co,e,o- 
The distribution is wider for (a) (0,0,0) = (0, ^,0) than for (b) (0,6,0) = (0, f ,0) and (c) (0,6,0) = (0, |§,0), showing the 
decrease in spread with increase in 6. The initial state of the particle is \^ins) = (|0) + i|l}) (8> IV'o) and the distribution is 
for 100 steps. 



corresponding to the quantum coin operation with £ = ( — and 9 — tt/A in Eq. {T}. The Hadamard operation is 
an unbiased coin operation, and the resulting walk is known as the Hadamard walk. This walk implemented on a 
particle initially in a symmetric superposition state, 

|* iras ) = -J=[|0)+i|l)]®|Vo), (7) 

obtained by choosing 8 — 7r/4 and r\ = ix jl in Eq. @, returns a symmetric, unbiased probability distribution of the 
particle in position space. However, the Hadamard walk on any asymmetric initial state of the particle results in an 
asymmetric, biased probability distribution of the particle in position space [3l|. We should note that the role of the 
initial state on the symmetry of the probability distribution is not vital for a QW using the three-parameter operator 
given by Eq. as a quantum coin operation. 

To elaborate this further, we will consider the first-step evolution of the DTQW using a three-parameter quantum 
coin operation given by Eq. Q on a particle initially in the symmetric superposition state. After the first step of the 
DTQW the state can be written as 

W SA i\*ins) = ^= [(e l «cos(#) +i e *sin(0)) |0)|V-i) + (e~* sin(0) - ie~^ cos(0)) |l)|V+i>] ■ (8) 

If £ = Q, Eq. ([5]) has left-right symmetry in the position probability distribution, but not otherwise. That is, the 
parameters £ and £ introduce asymmetry in the position space probability distribution. Therefore, a coin operation 
with £ ^ £ in Eq. ([T|) can be called as a biased quantum coin operation which will bias the QW probability distribution 
of the particle initially in a symmetric superposition state (Fig. [5J [3q . However, we should note that irrespective of 
the quantum coin operation used, QW can also be biased by choosing an asymmetric initial state of the particle (for 
example, the Hadamard walk of a particle initially in the state |0) or the state |1)). 

B. Many-particle quantum walk 

To define a many-particle QW in one dimension, we will consider an A/-particle system with one non-interacting 
particle at each position (Fig. [3]). The M identical particles in M lattice points with each particle having its own 
coin and position Hilbert space will have a total Hilbert space T-L = (H c ® T-L P ) M ■ We assume the particles to be 
distinguishable. 



-(a) = (0°, 15°, 0°) 
-(b) = (0°, 45°, 0°) 
(c) = (0°, 75°, 0°) 
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FIG. 2: Spread of probability distribution for different values of £, 0, £ using the quantum coin operator C^e,(' The parameter 
£ shifts the distribution to the left: (a)(£, 9, £) = (§ , f , 0) and (c) (£, 0, f) = |, 0). The parameter £ shifts it to the right: 
(b) (£,0,0 = (0, f , f ) and (d) (£,0,C) = (0, f , f§ ). The initial state of the particle |* lns ) = -^(|0) + ® |^o> and the 
distribution is for 100 steps. 



The evolution of each step of the QW on the M-particle system is given by the application of the operator W®q Q . 
The initial state that we will consider for the many-particle system in one dimension will be 

■_ M — l 

\0)+i\l) 



I0> shifts to left 
ll> shifts to right 



AAA A A A A A 





M particles and 2t+M+l lattice positions 
FIG. 3: Many-particle state with one non-interacting particle at each position space. 



For an M-particle system after t steps of the QW, the Hilbert space consists of the tensor product of single lattice 
position Hilbert space which is (2t + M + I) in number. That is, after t steps of the QW, the M particles are spread 
between (J — t) to (J + t). In principle, each lattice point is associated with a Hilbert space spanned by two subspaces, 
a zero-particle subspace and one-particle subspace spanned by two possible states of the coin, |0) and |1). Therefore, 
the dimension of each lattice point will be 3 M and the dimension of total Hilbert space is (3 M )® M . pjg Q shows 
the probability distribution of the many-particle system with an increase in number of steps of the QW. 



III. ENTANGLEMENT 



To efficiently make use of entanglement as a physical resource, the amount of entanglement in a given system has 
to be quantified. Therefore, entanglement in a pure bipartite system or a system with two Hilbert spaces is quantified 
using standard measures known as entropy of entanglement or Schmidt number The entropy of entanglement 
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FIG. 4: (color online) Probability distribution of 40 particles initially with one particle in each position space when subjected 
to the QW of different number of steps. The initial state of all the particles is -75(10} + and is evolved in position space 
using the Hadamard operator, as the quantum coin. The distribution spreads in the position space with an increase 

in number of steps. 



corresponds to the von Neumann entropy, a functional of the eigenvalues of the reduced density matrix, and a Schmidt 
number is the number of non-zero Schmidt coefficients in its Schmidt decomposition. For a multipartite state, there are 
quite a few good entanglement measures that have been proposed (37rl43| . However, as the number of particles in the 
system increases, the complexity of finding an appropriate entanglement measure also increases, making scalability 
impractical. Among the proposed measures, to address this scalability problem, Mayer and Wallach proposed a 
scalable global entanglement measure (polynomial measure) to quantify entanglement in many-particle systems [40j ] . 

In this section, we will first discuss the entanglement of a particle with position space quantified using entropy of 
entanglement. Later we will discuss spatial entanglement quantified using the Mayer- Wallach (M-W) measure. Spatial 
entanglement has been explored earlier using different methods. For example, in an ideal bosonic gas it has been 
studied using off-diagonal long-range order [44| . For our investigations, we consider a distinguishable many-particle 
system, implement QW and use the M-W measure to quantify spatial entanglement. In this system the dynamics of 
particles can be controlled by varying the quantum coin parameters, the initial state of the particles, the number of 
particles in the system and the number of steps of the QW. In particular, we choose the particles in one-dimensional 
open and closed chains. The spatial entanglement thus created can be used for example to create entanglement 
between distant atoms in an optical lattice [45| or as a channel for state transfer in spin chain systems [46l448l | . 



A. Single-particle - position entanglement 



QW entangles the particle (coin) and the position degrees of freedom. To quantify it, let us consider a DTQW on 
a particle initially in a state given by Eq. j7]) with a simple form of a coin operation 

_ / cos(0) sin(0) \ , , 

Co ^ 0= {sm(0) -cos(0) ) ■ (10) 

After the first step, Wo, 0,0 = S(C ^s.o ® 1), the state takes the form 

|*i) = Wb,fl,o|*« M >=7(|0)gi|Vi-i» + <y(|l>®l^+i» (11) 

where 7 = ^ cos ( 6 ')+y m ( 9 ) ^ and 5 = ^ Sln ( 9 )~y° s ( e ) ^ - The Schmidt rank of is 2 which implies entanglement in the 

system. The value of entanglement with an increase in the number of steps can be further quantified by computing 
the von Neumann entropy of the reduced density matrix of the position subspace. 
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FIG. 5: (color online) Entanglement of a single particle with position space when subjected to the QW. The initial state of 
a particle is -/g(|0) + and is evolved in position space using different values for 8 in the quantum coin operation Co,e,o- 
The entanglement initially oscillates and approaches an asymptotic value with an increase in the number of steps. For smaller 
values of 9 the entanglement is higher and decreases with an increase in 8. Initial oscillation is also larger for higher 8. 
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FIG. 6: (color online) Entanglement of single particle with position space when subjected to the QW. The initial state of the 
particle is given by Eq. ([2| with 5 = 22L and rj = ^ and is evolved in position space using different values for 8 in the quantum 
coin operation Co,e,o- The entanglement initially oscillates and approaches an asymptotic value with an increase in the number 
of steps. For smaller values of 8 the entanglement is higher and decreases with an increase in 8. Initial oscillation is also larger 
for higher 8. 



Fig. [5] shows a plot of the entanglement against the number of steps of the QW on a particle initially in a symmetric 
superposition state using different values for 9 in the operation We- The von Neumann entropy of the reduced density 
matrix of the coin is used to quantify entanglement between the coin and the position in Fig. [5j That is, 



E c (i) = -X>]og a (Aj) 



(12) 



where Xj are eigenvalues of the reduced density matrix of the coin after t steps (time). The entanglement initially 
oscillates and reaches an asymptotic value with increasing number of steps. In the asymptotic limit, the entanglement 



7 



value decreases with an increase in 9 and this dependence can be attributed to the spread of the amplitude distribution 
in position space. That is, with an increase in 9, constructive interference of quantum amplitudes toward the origin 
becomes prominent narrowing the distribution in the position space. In Fig. [6l the process is repeated for a particle 
initially in an asymmetric superposition state \^i n ) = [ cos (^f )|0) + e l % sin( 2 ^ L )|l}] ® |^>o)- Comparing Fig. [5] with 
Fig. [5j we can note the increase in entanglement and decrease in the oscillation. This observation can be explained by 
going back to our earlier note on biased QW in Sec. Ill Al In Fig. [2] we note that biasing of the coin operation leads 
to an asymmetry in the probability distribution, with an increase in peak height on one side and a decrease on the 
other side (increase and decrease are in reference to the symmetric distribution). A similar biasing effect can also be 
reproduced by choosing an asymmetric initial state of the particle. The biased distribution with an increased value of 
probability at one side in the distribution contributes to a reduced oscillation in the distribution. This in turn results 
in the increase of the von Neumann entropy: entanglement. 

B. Spatial entanglement 

Spatial entanglement is the entanglement between the lattice points. This entanglement takes the form of non-local 
particle number correlations between spatial modes. To observe spatial entanglement we first need to associate the 
lattice with the state of a particle. Then we need to consider the evolution of a single-particle QW followed by the 
evolution of a many-particle QW, in order to understand spatial entanglement. 

1. Using a single-particle quantum walk 

In a single-particle QW, each lattice point is associated with a Hilbcrt space spanned by two subspaces. The first 
is the zero-particle subspace which does not involve any coin (particle) states. The other is the one-particle subspace 
spanned by the two possible states of the coin, |0) and To obtain the spatial entanglement we will write the state 
of the particle in the form of the state of a lattice. Following from Eq. ([TTj) . the state of the particles after first two 
steps of QW takes the form 

|* 2 ) = Woao|*i>=7[cob(0)|O)|V^ (13) 

In order to obtain the state of the lattice we can redefine the position state in the following way: the occupied position 
state \tjjj} as |L,-), which means that the j th position is occupied and the rest of the lattice is empty. Therefore, we 
can rewrite Eq. (fT5|) as 

|* 2 ) = 7 [cos(0)|O)|l,_ 2 ) + S in(0)|l)|l i )] +*[sin(fl)|0)|l J -)-cos(fl)|l)|l i+2 )]. (14) 

Since we are interested in the spatial entanglement, we project this state into one of the coin state so that we can 
ignore the entanglement between the coin and the position state and consider only the lattice states. Here we will 
choose the coin state to be |0) and take projection to obtain the state of the lattice in the form 

|*i«t) = |0) (TcoB^lVa) + 5smW)\h)) • ( 15 ) 

Each lattice site j can be considered as a Hilbcrt space with basis states (occupied state) and |0j) (unoccupied 
state) . Then, the above Eq. ([15]) in the extended Hilbert space of each lattice can be rewritten in terms of occupied 
and unoccupied lattice states as 

|*{ at )=7cos(fl)|l i _ 2 J )+<ysiii(e)|0 J -_ 2 l J -). (16) 

We can see that after first two steps of the QW the lattice points j and (j — 2) are entangled. One can check that the 
lattice points j and (j + 2) are entangled if we choose the coin state to be |1). With an increase in the number of steps, 
the state of the particle spreads in position space and the projection over one of the coin state reduces that state to 
a pure state, for which one may compute spatial entanglement, according to the above prescription. Therefore, with 
an increase in the number of steps, the spatial entanglement from a single-particle QW decreases. 

2. Using many-particle quantum walk 



We will extend the study of evolution of spatial entanglement as the QW progresses on a many-particle system. 
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Let us first consider the analysis of first two steps of the Hadamard walk (6 — tt/4 in Eq. (JTOj) ) on a three-particle 
system with the initial state: 



+1 



i*£> = ® 



|0)+i|l) 



® IV»i>- 



(17) 



We will label the three particles at positions — 1, and 1 as A, B and C, respectively. Since evolution of these particles 
is independent, we write down the state after the first step as a tensor product of each of the three particles: 



w®3 |^3p v 



[7|0>|-2) + *|l)|0)] A ®[ 7 |' 



1 



5|1>| + 1)] B ®[7|0)|0)+5|1)| + 2)] C , 



(18) 



where 7 = (1 + i)/2 and S = (1 — i)/2. After two steps the tensor product of each of the three particles is given by 



/ |Q)|- 3) + |1)|-1) 

V V2 
|0>|-2) + |l)|0) 



+ 5 



-1) - ID 



V2 

1> + |1>I + 1> 



V2 

|0)|0)-|l)|+2) 



V2 
|0)| + 1) 



|l)l+3) 



V2 



(19) 



By projecting this state into one of the coin states (we choose state |0) <8> |0) <8> |0)) we can obtain a state of the lattice 
for which spatial entanglement may be computed. Then the state of the lattice after projection and normalization is 

|*i«t) =7 3 |A)- 3 |B)_ 2 |C)-i 



7 2 £( |A)- 3 |B)_ 2 |C)i + |A)_ 3 [B>o[C>_i + |AC)_i|B)_ 2 ) 
7 <5 2 ( |A)_ 3 |B)o|C)a + |AC)_i|B) + |A)_i|B)_ a |C)i) 
<5 3 |A)_ 1 |B) |C) 1 , 



(20) 



where A, B and C represent the particle labels and the subscripts represent the position labels. In a similar manner 
we can obtain |\&j tt t) for a system with a large number of particles. Then the next task is to calculate the spatial 
entanglement. 

IV. CALCULATING SPATIAL ENTANGLEMENT IN A MULTIPARTITE SYSTEM 

In a system with two particles, the state is separable if we can write it as a tensor product of individual particle 
states, and entangled if not. For a system with M > 2 particles, a state is said to be fully separable if it can be 
written as 



|V) = |<M ® ® • •• |<M 



(21) 



when k = M . \4>i) will then denote the state of the i th particle. When k < M a state is said to be partially entangled 
and when k = 1 the state will be fully entangled. 

Rather than using the von Neumann entropy to quantify multipartite entanglement of a given state p, one sometimes 
often prefers to consider purity, which corresponds (up to a constant) to linear entropy, that is the first-order term in 
the expansion of the von Neumann entropy around its maxima, given by 



E = 



1 [1 - 



(22) 



for a d-dimensional particle Hilbert space [49|. To quantify the entanglement of multipartite pure states, one measure 
commonly, used is the Meyer- Wallach (M-W) measure |4Q| . It is the entanglement measure of a single particle to the 
rest of the system, averaged over the whole of the system and is given by 



E 



MW 



d-l 



(23) 
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where L is the system size and pi is the reduced density matrix of the i th subsystem. The M-W measure does not 
diverge with increasing system size and is relatively easy to calculate. 

In a multipartite QW the dimension at each lattice point, after projection over one particular state of coin, is 2 M 
where M is the number of particles. Hence, the expression for entanglement will be 



E M w{\4>lat)) 



2 M -1 



1 - 



1 



t-H 



2t + M + 1 



(24) 



where t is the number of steps and pj is the reduced density matrix of j th lattice point. The reduced density matrix 
Pj can be written as 



(25) 



where \k) is one of the 2 M possible states available for a lattice point and pi can be calculated once we have the 
probability distribution of an individual particle on the lattice. 

Since we have M distinguishable particles, we have 2 M configurations depending upon whether a given particle is 
present in the lattice point or not after freezing the state of the particle. This set of configurations forms the basis for 
a single- lattice point Hilbcrt space. Now we can calculate pi, the probability of k th configuration of a particle in the 

j th lattice point as follows. Let us say aj is the probability of the i th particle to be or not to be in the j th lattice 
point depending on ij. If li is 1, then it gives us the probability of the particle to be in the lattice point. If U is 0, 

Hence, we can write 



then a,j is the probability of a particle not to be in the lattice point, that is, a*- 

P* 1 1 • 



1 



(26) 



Once we have the probability of each particle at a given lattice position, the spatial entanglement can be conveniently 
calculated. Since the QW is a controlled evolution, one can obtain a probability distribution of each particle over 
all lattice positions. In fact, one can easily control the probability distribution by varying quantum coin parameters 
during the QW process and hence the entanglement. 
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FIG. 7: (color online) Evolution of spatial entanglement with an increase in the number of steps of the QW for different number 
of particles in an open one-dimensional lattice chain. The entanglement first increases and with further increase in the number 
of steps, the number of lattice positions exceeds the number of particles in the system resulting in the decrease of the spatial 
entanglement. The distribution is obtained by implementing the QW on particles in the initial state -75GO} + i|l)) and the 
Hadamard operation Cq ,W4,o as quantum coin operation. 



Fig. [7] shows the phase diagram of the spatial entanglement using a many-particle QW. Data for the phase diagram 
were obtained numerically by subjecting the many-particle system with different number of particles to the QW with 
increasing number of steps. The quantity of spatial entanglement was computed using Eq. (|24[) . 
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FIG. 8: (color online) Quantity of spatial entanglement for 10 particles after 10 steps and 20 particles after 20 steps of the QW 
on a one-dimensional lattice using different values of 9 in the quantum coin operation Co,e,o- For (a) and (b), the distribution 
is for particles initially in the symmetric superposition state, — ^= ( [0) + and for (c) the particle's initial state is |0) (will be 
the same for state |1)). Quantity of entanglement is higher for 9 closer to and 7r/2 compared to the intermediate value. We 
note that the asymmetric probability distribution due to an asymmetric initial state in case of (c) contributes for an increase 
in the quantity of spatial entanglement. When 9 = 7r/2, for every even number of steps of the QW, the system returns to the 
initial state where entanglement is 0. Entanglement is for 9 — 0. 
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FIG. 9: (color online) Evolution of spatial entanglement for a system with different number of particles in a closed chain. With 
an increase in the number of steps, the entanglement value remains close to asymptotic value with some peaks in between. The 
peaks can be accounted for the crossover of leftward and rightward propagating amplitudes of the internal state of the particle 
during the QW. The peaks are more for a chain with a smaller number of particles. An increase in the number of particles in 
the system results in the decrease of the entanglement value. The distribution is obtained by using -^(|0) +i|l)) as the initial 
states of all particles and the Hadamard operation Cqw^o as quantum coin operation. 




— 10 particles 

— 14 particles 

— 20 particles 



Here, we have chosen the Hadamard operation Co i7r /4,o an d ^ (|0) + i|l)) as the quantum coin operation and initial 
state of the particles, respectively, for the evolution of the many-particle QW. To see the variation of entanglement 
for a fixed number of particles with an increase in steps, we can pick a line parallel to the y axis. That is, fix the 
number of particles and see the variation of entanglement with the number of steps. 

In Fig. we see that for a fixed number of particles, the entanglement at first increases to some value before 
gradually falling. For M = 12 we can note that the peak value is about 0.5 before gradually falling. With an increase 
in the number of steps of the QW, the number of lattice positions to which the particles evolve increases resulting in 
the decrease of the spatial entanglement (see Eq. (|24l) ). The decrease in entanglement before the number of steps is 
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FIG. 10: (color online) Quantity of spatial entanglement for 20 particles on a closed chain after 20 steps of the QW using 
different values of 8 in the quantum coin operation Co.e.o- The distribution is for particles initially in the state ^|(|0) + 
Since the system is a closed chain, the QW does not expand the position Hilbert space, and therefore for all values of 9 from 
to 7r/2 the entanglement value remains roughly uniform except for a small peak at smaller values of 9. For 9 — when the 
number of steps equal to the number of particles, the amplitudes goes round the chain and returns to its initial state making 
the entanglement and for 9 = tt/2, for every even number of steps of the QW, the system returns to the initial state where 
entanglement is again 0. 



equal to the number of particles should be noted. This is because for the Hadamard walk the spread of a probability 
distribution after t steps is between ^| and ^ (36j . 

If we fix the number of steps and measure the entanglement by increasing the number of particles in the system, 
the quantity of spatial entanglement first decreases and then it starts increasing with an increase in the number of 
particles. 

To show the variation of spatial entanglement with the quantum coin parameter 9, we plot the spatial entanglement 
by varying the parameter 9 for a system with 10 particles after 10 steps of the QW and for a system with 20 particles 
after 20 steps of the QW in Fig. [5] In this figure, (a) and (b) are plots that use the symmetric superposition state 
-^(|0) + (unbiased QW) as an initial state of all the particles, and (c) is the plot with all the particles in one 
of the basis states |0) or |1) (biased QW) as the initial state. We note that the quantity of entanglement is higher 
for 9 values closer to and ir/2 and dips for values in between for all the three cases. Biasing the QW, plot (c) 
shows a slight increase in the quantity of entanglement compared to the unbiased case, plot (b). A similar effect is 
seen by biasing the QW using two parameters £, £ in the coin operation C^g^ on particles initially in a symmetric 
superposition state. 

Closed chain: Since most physical systems considered for implementation will be of a definite dimension, we extend 
our calculations to one of the simplest examples of closed geometry, an n— cycle. For a QW on an n— cycle, the shift 
operation, Eq. ((4]), takes the form 

n — 1 n — 1 

S = |0)(0| ® Y, l^-i mod + ® E l^+i mod „>(^l- (27) 

3=0 3=0 

When we consider a many-particle system in a closed chain, with the number of lattice positions equal to the number 
of particles M, the QW process does not expand the position Hilbert space like it does on an open chain (line). 
Therefore the spatial entanglement does not decrease at later times as it does for a walk on an open chain, but 
remains close to the asymptotic value. Fig. [5] shows the evolution of entanglement for a system with different number 
of particles in a closed chain. The peaks seen in the plot can be accounted for by the crossover of the leftward and 
rightward propagating amplitudes of the internal state of the particle during the QW process. The frequency of the 
peaks is more for a smaller number of particles (smaller closed chain). Also, note that the increase in the number of 
particles and the number of lattice points in the closed cycle results in the decrease in spatial entanglement of the 
system. 
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In Fig. [TUl the value of spatial entanglement for 20 particles on a closed chain after 20 steps of the QW using different 
values of 9 in the quantum coin operation Co,e o is presented. For all values of 9 from to it/2 the entanglement 
value remains roughly uniform except for the extreme values of 9. For = 0, the amplitude goes round the ring and 
returns to its initial state making the spatial entanglement value = 0. For 9 = 7r/2, for every even number of steps of 
the QW, the system returns to the initial state where spatial entanglement is again 0. 

Therefore, spatial entanglement on a large lattice space can be created, controlled and optimized for a maximum 
entanglement value by varying the quantum coin parameters and number of particles in the multi-particle QW. 



V. CONCLUSION 



We have presented the evolution of spatial entanglement in a many-particles system subjected to a QW process. 
By considering many particle in the one-dimensional open and closed chain we have shown that spatial entanglement 
can be generated and controlled by varying the quantum coin parameters, the initial state and the number of steps in 
the dynamics of the QW process. The spatial entanglement generated can have a potential application in quantum 
information theory and other physical processes. 
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